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Abstract 

We establish sublinear growth of correctors in the context of stochas¬ 
tic homogenization of linear elliptic PDEs. In case of weak decorre¬ 
lation and “essentially Gaussian” coefficient fields, we obtain optimal 
(stretched exponential) stochastic moments for the minimal radius 
above which the corrector is sublinear. Our estimates also capture 
the quantitative sublinearity of the corrector (caused by the quanti¬ 
tative decorrelation on larger scales) correctly. The result is based on 
estimates on the Malliavin derivative for certain functionals which are 
basically averages of the gradient of the corrector, on concentration of 
measure, and on a mean value property for a-harmonic functions. 


1 Introduction 

In the present work, we are concerned with the stochastic homogenization of 
linear elliptic equations of the form 


-V • a\/u = f. (1) 

In stochastic homogenization of elliptic PDEs, a is typically a uniformly 
elliptic and bounded coefficient field, chosen at random according to some 
stationary and ergodic ensemble (•). On large scales (and for slowly varying 
/), one may then approximate the solution u to the equation (JTJ) by the 
solution Uhom to the so-called effective equation 


V ' d/iomV Uh om f , 


( 2 ) 
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which is a constant-coefficient equation with the so-called effective coeffi¬ 
cient ahom- Mathematically, this homogenization effect is encoded in growth 
properties of the corrector (cf. below for a definition of the corrector). 

The goal of the present paper is to provide a fairly simple proof of quanti¬ 
fied sublinear growth of the corrector under very mild assumptions on the 
decorrelation of the coefficient held a under the ensemble (•). We do this in 
the context of coefficient fields that are essentially Gaussian. More precisely, 
we consider coefficient fields a which are obtained from a Gaussian random 
held by pointwise application of some (nonlinear) mapping, the role of the 
nonlinear map being basically to enforce uniform ellipticity and boundedness 
of our coefficient held. 

The motivation for this result is the following: Gloria, Neukamm, and the 
second author ra have shown that qualitatively sublinear growth of the 
(extended) corrector (</>, a) (cf. below for a definition) entails a large-scale in¬ 
trinsic C 1,a regularity theory for a-harmonic functions. In a subsequent work 
[8], the two authors of the present paper have shown that slightly quantified 
sublinear growth of the corrector even leads to a large-scale intrinsic C k,a 
regularity theory for any k G N. Therefore, the results of the present work 
show that even in case of ensembles with very mild decorrelation, for almost 
every realization of the coefficient held, a-harmonic functions have arbitrary 
intrinsic smoothness properties on large scales. Furthermore, our results en¬ 
able us to estimate the scale above which this happens - a random quantity - 
in a stochastically optimal way. Indeed, the motivation for the present work 
was to establish such a (necessarily intrinsic) higher order regularity theory 
under the weakest possible assumptions on the decay of correlations. 

By an intrinsic regularity theory we mean that the regularity is measured 
in terms of objects intrinsic to the Riemannian geometry defined by the 
coefficient held a, like the dimension of the space of a-harmonic functions of 
a certain algebraic growth rate, or like estimates on the Holder modulus of the 
derivative of a-harmonic functions as measured in terms of their distance to 
a-linear functions. An extrinsic large-scale regularity theory for a-harmonic 
functions in case of random coefficients was initiated on the level of a C 0,a 
in [5J [18] and pushed to C 1,0 in [3], which significantly extended qualitative 
arguments from the periodic case [1] to quantitative arguments in the random 
case. However, an extrinsic regularity theory is limited to C 1,0 , as can be 
seen considering the harmonic coordinates: Taking higher order polynomials 
into account does not increase the local approximation order. 

After this motivation, we now return to the discussion of the history on 
bounds on the corrector, as they depend on assumptions on the stationary 
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ensemble of coefficient fields. Almost-sure sublinearity (always meant in a 
spatially averaged sense) of the corrector 0 under the mere assumption of 
ergodicity was a key ingredient in the original work on stochastic homoge¬ 
nization by Kozlov [16] and by Papanicolaou & Varadhan m- Almost-sure 
sublinearity of the extended corrector (0, er), as is needed for the large-scale 
intrinsic C 1 ’"-regularity theory, was established in ra under mere ergodicity. 

Yurinskii j22| was the first to quantify sublincar growth under general mix¬ 
ing conditions, however only capturing suboptimal rates even in case of finite 
range of dependence. Very recently, a much improved quantification of sub- 
linear growth of 0 under finite range assumptions was put forward by Arm¬ 
strong, Kuusi, and Mourrat p], relying on a variational approach to quanti¬ 
tative stochastic homogenization introduced by Armstrong and Smart [3] , an 
approach which presumably can be extended to the case of non-symmetric 
coefficients and more general mixing conditions following [2], However, while 
this approach gives the optimal, i. e. Gaussian, stochastic integrability, it 
presently fails to give the optimal growth rates. 

Optimal growth rates have been obtained under a quantification of ergodic¬ 
ity different from finite range or mixing conditions, namely under Spectral 
Gap assumptions on the ensemble. This functional analytic tool from sta¬ 
tistical mechanics was introduced into the field of stochastic homogenization 
in an unpublished paper by Naddaf and Spencer [20], and further leveraged 
by Conlon et. ah yielding optimal rates for some errors in stochastic 

homogenization in case of a small ellipticity contrast. The work of Gloria, 
Neukamm and the second author extended these results to the present case 
of arbitrary ellipticity contrast mmm, in particular yielding at most log¬ 
arithmic growth of the corrector (and its stationarity in d > 2). Loosely 
speaking, the assumption of a Spectral Gap Inequality amounts to correla¬ 
tions with integrable tails; in the above-mentioned works it has been used 
for discrete media (i. e. random conductance models), but has subsequently 
been extended to the continuum case mm- 

A strengthening of the Spectral Gap Inequality is given by the Logarith¬ 
mic Sobolev Inequality (LSI); it is a slight strengthening in terms of the 
assumption (still essentially encoding integrable tails of the correlations), 
but a substantial improvement in its effect, since it implies Gaussian concen¬ 
tration of measure for Lipschitz random variables. The assumption of LSI 
and implicitly concentration of measure, which will be explicitly used in this 
work, has been introduced into stochastic homogenization by Marahrens and 
the second author [18] . In ra. it has been shown that the concept of LSI can 
be adapted to also capture ensembles with slowly decaying correlations, i. e. 
thick non-integrable tails, by adapting the norm of the vertical or Malliavin 
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derivative to the correlation structure. As a result, the stochastic integra- 
bility of the optimal rates could be improved from algebraic to (stretched) 
exponential, but missing the expected Gaussian integrability. 

The main merit of the present contribution w. r. t. to pjOj is twofold: First, our 
approach directly provides optimal quantitative sublinearity of the corrector 
(0, a) on all scales above a random minimal radius r*, i.e. in contrast to 
the estimates of ra our estimates capture the decorrelation on scales larger 
than r* in a single argument. Note that our definition of r* differs from the 
one in [10J. Second, in case of weak decorrelation, our simpler arguments 
are nevertheless sufficient to establish optimal stochastic moments for the 
minimal radius r* above which the corrector (0, cr) displays the quantified 
sublinear growth. 

In the present work, we consider the following type of ensembles on A- 
uniformly elliptic tensor fields a = a(x ) on M. d : Let a = a(x) be a tensor¬ 
valued Gaussian random held on W l that is centered (i. e. of vanishing ex¬ 
pectation) and stationary (i. e. invariant under translation) and thus charac¬ 
terized by the covariance (a(x) <S> a(0)). Our only additional assumption on 
a is that there exists an exponent f3 G (0, d) such that 

|(a(a;) <g) a(0))| < |a:| _/3 for all x G M. d . (3) 

In this work, we are concerned with the case of weak decay of correlation 
in the sense of /3 -C 1. Let $ be a 1-Lipschitz map from the space of 
tensors into the space of A-clliptic symmetric tensors. Then our ensemble 
is the distribution of a where a is given by a(x) := <h(a(a;)). Note that the 
normalization in the constant in (J3]) and in the Lipschitz constant is not 
essential, since it can be achieved by a rescaling of x and the amplitude of a. 

Concerning the mathematical tools of our approach, several ideas are inspired 
by the work HUj. In particular, a key component of our approach are sensitiv¬ 
ity estimates (Malliavin derivative bounds) for certain integral functionals, 
which basically average the gradient V(0, cr) over an appropriate cube. Fur¬ 
thermore, we rely on a mean-value property for a-harmonic functions, which 
has been derived in m under appropriate smallness assumptions on the 
corrector. In our present contribution, we however pursue a conceptually 
simpler route to estimate the Malliavin derivative: The sensitivity estimate 
is performed through appropriate lA-norm bounds and Meyer’s estimate, 
rather than a more involved £ 2 — Id-norm bound like in |10j . 

Before stating our main results, let us recall the concept of correctors in 
homogenization and introduce some notation. The basic idea underlying 
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the concept of correctors in homogenization is the observation that the os¬ 
cillations in the gradient of solutions to the homogenized (constant- 

coefficient) problem (J2]) occur on a much larger scale than the oscillations 
in the gradient Vw of solutions to the original problem (ITj) . Thus, it is im¬ 
portant to understand how to add oscillations to an affine map (an affine 
map being always n/,, OTO -harmonic) to obtain an a-harmonic map. In the con¬ 
text of stochastic homogenization, one is therefore interested in constructing 
random scalar fields (pi = (pi(a,x ) subject to the equation 

- V ■ a(ei + X7(j)i ) = 0 (4) 

which almost surely display sublinear growth in x: The (pi then facilitate 
the transition from the ahom~ harmonic (Euclidean) coordinates x Xi to 
the “a-harmonic coordinates” x K > Xi + (f>i(x). Since any affine map may be 
represented in the form b + JT p t x t for b , G M, the (pi also facilitate the 
construction of associated a-harmonic “corrected affine maps” b + ^2 i ^ i (x i + 
(pi)- 

With the help of the corrector, one may characterize the effective coefficient 
ahom'- In our setting of stochastic homogenization, the effective coefficient is 
given by the formula 

Q’hom&i {®(Ci T ^(pi))i (5) 

where (•) refers to the expectation with respect to our ensemble (i.e. proba¬ 
bility measure). 

In the language of a conducting medium with conductivity tensor a - note 
that in this picture, one has / = 0 in (pQ) - , the quantity E) := e* + Vd* 
corresponds to the (curl-free) “microscopic” electric field associated with a 
“macroscopic” electric field e* (and, therefore, (pi corresponds to the “mi¬ 
croscopic” correction to the “macroscopic” electric potential xp). The corre¬ 
sponding (divergence-free) “microscopic” current density is given by 

ft \= a(e-i + V&), (6) 

while the “macroscopic” current density associated with the “macroscopic” 
electric field e* is given by the “average” of this quantity, i.e. by the expression 

©• 

In periodic homogenization of linear elliptic PDEs, it turns out to be con¬ 
venient to introduce a dual quantity to the corrector (pi (cf. e.g. [T5j p.27]): 
One constructs a tensor field crskew-symmetric in the last two indices, 
which is a potential for the flux correction — ah om ei in the sense 

V • <Tj a(ei T V (pi) (7) 
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where we have set (V • cr,)^ := Ylk=i dk&ijk- With the help of this “extended 
corrector” (0, a), it is possible to give a bound on the homogenization error 
(in terms of appropriate norms of 0 and a). 

One of the main merits of [10] is the discovery of the usefulness of this 
extended corrector (0, a) in the context of stochastic homogenization. For 
stationary and ergodic ensembles (•) of A-uniformly elliptic and symmetric 
coefficient fields a = a(x) on K d , in phj correctors 0* and (Jijk such that 

are stationary, 

V0j,Voy / 'fc of bounded second moment, (8) 

and of vanishing expectation, 

have been constructed. As a consequence of this and of ergodicity, the 0j 
and (Tjjfe almost surely display sublinear growth. Note that in case of cq, the 
choice of the appropriate gauge is important for the property (jHJ) and for our 
work, as the equation (J7J) determines cq (which by its skew-symmetry and its 
behavior under change of coordinates may be identified with ad — 1-form) 
only up to the exterior derivative of a d — 2-form. In fact, the choice of the 
gauge in in is such that 

A (J'jj k — 9jC[ik dkQiji (9) 

which in view of (j4j) and (J6]) is clearly compatible with (171) . 

Notation. To quantify the cllipticity and boundedness of our coefficient 
fields, throughout the paper we shall work with the assumptions 

av ■ v > A|u| 2 for all v e R d , (10) 

|gw| < |v| for all v G M d , (11) 

where A G (0,1). Note that in view of rescaling, the upper bound (1111) on a 
does not induce a loss of generality of our results. 

For our convenience, throughout the paper we shall assume that our coeffi¬ 
cient held a is symmetric. The arguments however easily carry over to the 
case of non-symmetric coefficient fields by simultaneously considering the 
correctors for the dual equation (i.e. the PDE with coefficient held a*, a* 
denoting the transpose of a). 

The expression s < f is an abbreviation for s < Ct with C a generic constant 
only depending on the dimension d, the exponent 0 > 0, and the ellipticity 
ratio A > 0. 

The expression s<( stands for s < with C a generic sufficiently large 
constant only depending on the dimension d, the exponent f3 > 0, and the 
cllipticity ratio A > 0. 
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By 1(E) we denote the characteristic function of an event E. 

The notation f A f refers to the average integral over the set A, i.e. we have 

j~A f dx = §A f dx / L 1 dx - 

In the sequel, (0, cr) stands for any component 0*, (Tijk for i,j,k = 1, • • • , d. 

2 Main Results and Structure of Proof 

Let us now state our main theorem. To quantify the sublinear growth of 
the extended corrector (0, cr), we first quantify the decay of spatial averages 
of V(0, a) over larger scales. In view of the decorrelation assumption ([3]) 
for our ensemble of coefficient fields, we expect that, up to logarithms, it 
is the exponent | that governs the decay of averages of V(0, cr) and the 
improvement over linear growth for (0, cr). Indeed, this exponent is reflected 
in the theorem. 

Theorem 1. Let a = a(x) be a tensor-valued, Gaussian random field on 
that is centered (i. e. of vanishing expectation) and stationary (i. e. invariant 
under translation); assume that the covariance of a satisfies the estimate 

|(a(x) <g) a(0))| < |x| _/3 forallxeW 1 (12) 

for some fi € (0, d ). Let $ : W lxd —> M dxd be a Lipschitz map with Lipschitz 
constant < 1; suppose that <f> takes values in the set of symmetric matrices 
subject to the ellipticity and boundedness assumptions (fTU|h (1TTT) . Define the 
ensemble (•) as the probability distribution of a, where a is the image of a 
under pointwise application of the map <3>, i.e. a(x) := $(a(a;)). 

Assume in addition on the ensemble (•) that fi in (ED is sufficiently small in 
the sense of 

P<f, ( 13 ) 

where C denotes a generic constant only depending on d and X. 

i) Consider a linear functional F = Fh on vector fields h = h(x) satisfy¬ 
ing the boundedness property 

. r 2d . 

\Fh\<(4 \h\—pdx) 2d (14) 

J |#|<r 

for some radius r > 0. Then the random variable FV(0, a) satisfies 
uniform Gaussian bounds in the sense of 

(/(|FV(0, cr)| > M)) < C'exp(-^r /3 M 2 ) for all M < 1. (15) 

o 
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ii) There exists a (random) radius r* for which the “iterated logarithmic” 


bound 



forr > r * 
(16) 


holds and which satisfies the stretched exponential bound 


(exp(irf)) < C. 


( 17 ) 


Morally speaking, Theorem [T] converts statistical information on the coef- 
hcient held a (or rather a) into statistical information on the coefficient 
held V0 := V(0i,--- , fifi) related by (0J) - Despite the nonlinearity of the 
map a i —> V</>, which only in its linearization around a = id turns into 
the Helmholtz projection, Theorem Q] states that V</> essentially inherits the 
statistics of a: (1151) implies in particular that spatial averages F = jj x , <r Vfidx 
of V0 satisfy the same bounds as if V</> itself was Gaussian with correlation 
decay ([3]) . On the level of these Gaussian bounds, the only prize to pay for 
nonlinearity is the restriction M < 1 in (fT5l) on the threshold. 

Incidentally, the way we obtain ii) from i) bears similarities with an argument 
in [IJ in the sense that a decomposition into Haar wavelets is implicitly used. 

To obtain an estimate like (fT5l) . the starting point of our proof is the Gaussian 
concentration of measure applied to a. 

Proposition 1 (Concentration of Measure, cf. e.g. |17J Proposition 2.18])). 
Let a = a(x) be a tensor-valued Gaussian random field on M. d that is centered 
and stationary; denote its covariance operator by Gov. Consider a random 
variable F, that is, a function(al) F = F{a ). Suppose that F is 1-Lipschitz in 
the sense that its functional derivative, or rather its Frechet derivative with 
respect to L 2 (M d ; R dxd ), ^ = ^-{a, x), which can be considered a random 
tensor field and assimilated with a Malliavin derivative, satisfies 



for almost every a. (18) 


Then F has Gaussian moments in the sense of 


(exp (M(F — (F)))) < exp(-^-) for all M > 0. 


(19) 


Furthermore, for any M > 0 we have the estimate 
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( 20 ) 



We now substitute our assumption (fTBD on the Frechet derivative by a stronger 
but more tractable condition. 


Lemma 1. Let a = a(x ) be a tensor-valued, Gaussian random field on W 1 
that is centered and stationary; denote its covariance operator as Cov and 
suppose that for some /3 G (0, d ) we have the bound 


| (a(x) <g) a(0))| < |x| 13 forallx&W d . 

Let <f> : M. dxd —> M dxd be a 1-Lipschitz map; denote the probability distribution 
of <F(a) as (•). Consider a functional F on the space of tensor fields a of the 
form F = F(a) with a(x) := <F(a(a;)); we shall use the abbreviation F{a ) for 
F(<F(a)). Let q G (1,2) be given by 


1 

q 



( 21 ) 


and suppose that the Frechet derivative of F with respect to L 2 (M d ; M dxd ) 
satisfies 

^ J\-g—\ q dcc S j q 1 for (-)-almost every a. (22) 

Then the estimate m is satisfied, i.e. we have 



(a, x)(Cov 


dF 

da 


(a, -))(x)dx < 


1 


for almost every a. 


We observe that if q and (5 are related by (l2T]h as fi t d we have q f 2 and 
for fi J, 0 we have q f 1. 

For linear functionals of (the gradient of) the corrector (which are therefore 
nonlinear functionals of the coefficient held a), we now establish an explicit 
representation of the Frechet derivative; this will aid us in verifying the Lip- 
schitz condition (l22j) and thus ultimately the concentration of measure state¬ 
ments (1T9T) and (1201) for (an appropriate modification of) such functionals. 

p , , 

Lemma 2. Consider a linear functional on Lp zrT (M“;M“) of the form 


Fh := 


g ■ hdx , 


(23) 


where g G L p (R d ; M d ), p >2, and supp g C {|:c| < r} for some r > 1. Let a 
be some coefficient field subject to the ellipticity and boundedness conditions 
m, m- Then the following two assertions hold: 
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1) Consider the Frechet derivative ^ of the functional F := FV a, ]k (note 
that this functional is nonlinear in a, although it is linear in a. ljk ) at 
a (for some fixed i,j, k). Introduce the decaying solutions v, Vj k to the 
equations 

- Av = V • g (24) 

and (where a* denotes the transpose of a) 

- V • a*(Vvj k + (djVe k - d k vef)) = 0. (25) 

We then have the representation 
dF 

— (a) = ( djve k - d k vej + Vv jk ) ® (V<& + ef). (26) 

2) Consider the Frechet derivative 'jj( of the functional F := FVfi at a. 
Introduce the decaying solution v to the equation (again, a* denoting 
the transpose of a) 

- V • a*X7v = V • g. (27) 

We then have the representation 

dF 

—-(a) = Wv ® (V<f>i + ef). (28) 

oa 

The previous explicit representation of the Frechet derivative for certain lin¬ 
ear functionals of (the gradient of) the corrector (0, a) enables us to verify 
the bound (1221) for the Malliavin derivative, provided that a certain mean 
value property is satisfied for a-harmonic functions. Note that the latter 
requirement is a condition on the coefficient field a; in Lemma [4] below we 
shall provide a sufficient condition for this property to hold. 

As the functionals which the next lemma shall be applied to are basically 
averages of V0 or Va over cubes of a certain scale r, we state the lemma in a 
form which makes it directly applicable in such a setting. In particular, the 
boundedness assumption (J2H]) for the linear functional is motivated by these 
considerations. 

Lemma 3. Let r > 1 and consider a linear functional h i— y Fh on L~ (My M d ) 
(with p e (2, oo )) satisfying the support and boundedness condition 

\Fh\<(l \h\i^dx) 2 ^. (29) 

J |x|<r 
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Suppose that the constraint 


2 < p < 2 + c(d, A) 


(30) 


holds (with c(d, A) > 0 to be fixed in the proof below). Let q G (1, 2) be related 
to p through 



(31) 


Consider the Frechet derivative of the functional F := FX7<7ijk (or the 
functional F := FV0j/ note that these functionals are nonlinear functionals 
of a) at some symmetric coefficient field a subject to the conditions m, 

am ■ 

Provided that the coefficient field a is such that the mean value property 


j- \\7u\ 2 dx < + \\7u\ 2 dx for any R>r and any p e [r,R] (32) 

J\ X \<P J\x\<R 

holds for any a-harmonic function u and provided that furthermore a is such 
that 

.j iiri 4 (/ 1(0,^)-/ 

R-^oo R \J\ X \< R J | x | 

is satisfied, we have the estimate 

/ f 8F \ f (p-2)d 

(J \-gff\ qdx ) ~ r P ■ ( 34 ) 

Note that for q related to (3 through (Ell) and p related to q through (Eli) . 

_ (p — 2)d ' 

we have r p = r _|0 , i.e. by (I34p the L q - norm of the Malliavin derivative 
decays like r~%. This demonstrates that for functionals like our averages 
of V (0, a) - note that these functionals have vanishing expectation due to 
the vanishing expectation of V (0, a) - , the concentration of measure indeed 
improves on large scales with the desired exponent: The “typical value” of 
the average of V(0, a) on some scale r decays like . 

We now have to provide a sufficient condition for the mean value property for 
a-harmonic functions ([32]) . To do so, we make use of the following result from 
pm, which provides the mean-value property assuming just an appropriate 
sublinearity condition on the corrector (0, cr). 

Proposition 2 (see m Lemma 2]). There exists a constant Co only depend¬ 
ing on dimension d and ellipticity ratio A > 0 with the following property: 


|<i? 


(0W)| : 


X 1/2 

dx j =0 


(33) 
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Suppose that for an elliptic coefficient field a subject to the ellipticity and 
boundedness conditions m and m the scalar and vector potentials (0, a), 
cf. and satisfy 


(/ |(0,a)-/ (0, cr)\ 2 dx) ^ 

x J\x\<R J\x\<R 7 

Then for any two radii R> r and p e [r, R] 
in {|x| < R} we have 


< pg-R for all R> r. (35) 
and any a-harmonic function u 



\Vu\ 2 dx < 



\Wu\ 2 dx. 


We shall show in the proof of the next lemma that the quantitative sublinear- 
ity condition on the corrector (l35]i may be reduced to a smallness assumption 
on a certain family of linear functionals of the gradient of the corrector. Ba¬ 
sically, these functionals will be obtained by averaging the gradient of the 
corrector on appropriate cubes, cf. the proof below. In combination with the 
previous proposition, we get the following lemma. 

Note that this result will allow us to buckle, since by Lemma [3] and Lemma Q] 
the mean-value property (1321) and thus ultimately (T36|) implies 


,dF, 


n,R 


da 


< R-P. 


Lemma 4. There exist N < 1 linear functionals iq, ■ • - , F N satisfying the 
support and boundedness condition 

\Fnh\<{[ \h\^dx) P -^ 

J\x\<l 

such that the following holds: Denote by F n) R the rescaling of F n given by 
F n jfii = Let r > 0. Provided that the condition 

sup -Fn,i?V(0, a) < 1 (36) 

R>r dyadic\n=lp” ,7V 


and the condition (l33|i are satisfied, we have the smallness estimate (ESD for 
the corrector; in particular, by Proposition d the mean value property (EI|) 
holds for any a-harmonic function u on scales > r, i.e. we have 



\Vu\ 2 dx < / \Vu\ 2 dx 
J\x\<R 


for any R> r and any p G [r, R]. 
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With these preparations, we are able to establish our main theorem. The 
main technical difficulty in the proof below is that our estimate 



for the Malliavin derivative of linear functionals of the gradient of the cor¬ 
rector (cf. (EU) is a conditional bound: It relies on the assumption that 
the mean-value property (l32]i holds for a-harmonic functions on scales larger 
than r. For the concentration of measure estimate (USD, however, an uncon¬ 
ditional estimate of the form (TTBP or (J22D (the latter being a proxy for ([TSD ) is 
needed. By Lemma H] we know that the mean-value property holds, provided 
that for a certain family of linear functionals of the corrector the smallness 
estimate 

sup F ntR V((j),cr) < ^ r 

R>r dyadic;n=l,-" ,N V'O 

is satisfied (Co being a universal constant). To circumvent this problem, in 
the proof below we therefore introduce the family of functionals 

F r \= min { sup (0, a ), -J- ), 

^ R>r dyadic;n= I,-** ,7V ^0 ' 

for which by design the unconditional bound for the Malliavin derivative 



holds. Therefore, concentration of measure is applicable to F r . The remain¬ 
der of the proof of the first part of our theorem below is dedicated to handling 
the (a priori unknown) expectation (F r ). 

The proof of the second assertion of our main theorem will mainly rely on 
the first assertion of the theorem as well as the quantitative improvement of 
the Malliavin derivative of averages of (V0, Ver) on larger scales, as captured 
by the estimate (l34lh 

3 Concentration of Measure and Estimates of 
the Malliavin Derivative 

3.1 Concentration of Measure 

Proof of Proposition Ql For the proof of the concentration of measure esti¬ 
mate (USD, we refer the reader to [13 Proposition 2.18]. We now estab¬ 
lish (l2np . By Chebychev’s inequality, implies ( I(F — (F ) > M )) < 


13 



exp(—^§-). In combination with the same estimate with F replaced by —F, 
we obtain flUD- □ 

Proof of Lemma [3 We need to verify that the condition (ITS]) is implied by 
the assumption (I22]h 

To do so, we first note that by Holder’s inequality we have for any exponent 
1 < q < oo 


dF dF 
—Cov—cte < 
da da 



in 9F 

|Cov aa 



9-1 

9 


Since Cov is the convolution with (a(x)a(O)) and since we have the bound 
|(a(x) <g) a(0))| < \x\~P, we have for the second factor 


| Cov 


OF , \^ 


da 


! ‘~ x dx 


< 


3F 


\x — y\P da 


(: v)\dy 
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9-1 


dx 


■ 7—1 

9 


which allows us to use the Hardy-Littlewood-Sobolev inequality 


x 


, dF 


y\P da 


(■ y)\dy 


9 

9-1 


dx 


9-1 


< 


dF 

da 


\ q dx)\ 


provided the exponents q and f3 are related by f!2ip . From this string of 
inequalities we learn that (TT8|) also holds provided 


dF, 

da 


| q dx) -Cl for almost every a. 


(37) 


We now change variables according to a(x) = d>(a(a;)); by the chain rule 
for F(a) = F(d>(a)) we have a,x ) = $'(a(x))^(a, x), so that by the 1- 
Lipschitz continuity of <F, our assumption UTZb implies (1571) and thus (ITS]) . □ 


3.2 Representation of the Malliavin derivative 

Proof of Lemma 0 We first give the argument for the “vector potential” cr, 
fixing a component a t jk- Consider a functional of the form F := FWa^k 
with Fh as in (]23|) . We claim that the Frechet derivative of F with respect 
to a is given by (126|) where the functions v = v(x) and Vjk = Vjk(a,x ) are 
determined as the decaying solutions of the elliptic equations (124|) and (]25l) . 
Computing the functional derivative of F as a function of a amounts to a 
linearization. We thus consider an arbitrary tensor field 8a = 8a(x), which 
we think of as an infinitesimal perturbation of a, and which thus generates 
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infinitesimal perturbations 8<f and 5 a of f> and a according to (HD, and 

flSD, that is, 

— V ■ ( aV5(j>i + 8a(V(j)i + ef)) = 0 (38) 

and 

- A 8a ijk = dj[8a(W(j)i + ef) + aV(5^) fc - <9 fc (<5a(V0j + ef) + aVSfa).. (39) 

In terms of the infinitesimal perturbation 8F of F, this implies by integration 
by parts (or rather by directly appealing to the weak Lax-Milgram formula¬ 
tions of the elliptic equations) 


8F = 


g ■ V Saijkdx 


I Vv ■ V 8aijkdx 

(djVek — dkvej) ■ (<5a(V</>j + ef) + aVSfifjdx 
(djvek — dkvej ) • <5a(V0j + ef)dx — J V'djk ■ a'V8(j) i dx 
(djVek — d^ej + Vhjfc) ■ 8a(V(j)i + efjdx, 


which is nothing else than (f26l) . 

Let us now establish the second part of our lemma. Consider a functional 
of the scalar potential of the form F := F'Vfa. To represent its Frechet 
derivative, introduce the decaying solution v to the equation (j2Tft . We observe 
that the variation of F with respect to a is given by 


8F = J g ■ VScpidx 
?P - [ a *X7v ■ VS&dx 


\7v ■ 8a(V(j)i + ei)dx , 


which leads to the conclusion 


□ 


3.3 Sensitivity estimate 

Proof of Lemma 0 We now argue that under certain boundedness assump¬ 
tions on F = Fh as a linear functional in vector fields h = h(x), we control 


15 


the size (1U2D of its Frechet derivative QF = Qj^(a,x) as a nonlinear func¬ 
tional FVcjjjfc = F(a') in coefficient fields a — a(x) (and similarly in the case 
F(a) = FV(f>i) for this case, the (simpler) proof is sketched afterwards). 

To this aim, let us first note that we have a Calderon-Zygmund estimate for 
—V • ftV with the exponents p and its dual exponent For any decaying 
function w and vector held h on W l related by 


— V ■ aVrc — V-h 


(40) 


we have 





and 


I\X7w\ p dx < 


\h\ p dx. 


(41) 


This assertion holds by Meyer’s estimate (see e.g. [19]), which only requires 
the cllipticity and boundedness assumptions (ITU]) . (ITT]) on a as well as the 
estimate \p — 2| <C 1, which is ensured by our condition (I3U]) . Note that an 
analogous estimate would hold for the dual equation —V • a*Vw = V ■ h .if 
our coefficient held were nonsymmetric. 

In the following, we will use the abbreviation || • || Pi # for the spatial L p -norm 
on the set B ; we write || • || p when B = K°\ We start by arguing that because 
G (1, 2), (l32|) also entails 





p- x dx. 


(42) 


It is obviously enough to establish (|42|) only for R > 2p; hence by Jensen’s 
inequality, (1T21) follows from (1321) once we establish the reverse Holder in¬ 
equality 


(4 \Vu\ 2 dx) 


1 

2 < 
r^j 


\Vu\dx. 


(43) 


J\x\<R 

To this purpose, we test —V • a\7u = 0 with r] 2l {u — m), where p is a smooth 
cut-oh of in {|x| < R} (with the property |V?y| < -^) and where 

the exponent 7 > 1 and the constant m G M will be chosen later. By the 
ellipticity and boundedness assumptions (fTU|) . (ITT]) and Young’s inequality 
we obtain 

[ (? 7 7 |Vu|) 2 cte < [((u — m)\'Vp' y \) 2 dx, 


and thus 


|V( 7 7 (m — m))\ 2 dx < / ((u — m)\X7if\) 2 dx, 
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which by the estimate on Vr/ gives 


||V(r 7 7 (w-m ))|| 2 < ^\\rp 1 (u-m) || 2 . (44) 

On the r. h. s. of (PHD we use first Holder’s inequality, then the isoperimetric 
inequality on {|ccr| < R} and finally Sobolev’s inequality on the whole space 
(for simplicity, we assume d > 2 here) 

\W~\u-m)\\ 2 < Wiu-m^l^Wu-mW^ 

d -2 d-1 ’l x l— 

< l|V(V(«-m))|| 2 1?1 ||V«.||’ |l|SR , 

provided the exponent 7 G (1, 00 ) is chosen such that \ — (1 — ^)^p + 

(which - as a simple computation shows - is satisfied precisely for 7 = |) and 
the constant m is the spatial average of u on {|x| < i?}. The combination of 
the last two estimates yields 

||V(V(« - m))|| 2 < 4||V(V(« - m))\\F ||V«||’ WSS , 

which (by 7 = f) entails ||Vw|| 2iN <| < ||V(rf (w - m ))|| 2 < Rri ||'Vm||i, n < r 
and thus (l43|) . 

We now give the argument for (134|) in case of a functional of the form 
(the case FVcfii will be treated below). Clearly fl29]i implies that there exists 
a (deterministic) vector field g — g(x) with 

suppg C {|cc| < r} and \\g\\ p < r~~ (45) 

such that we have the representation for F = Fh as a linear functional on 
vector fields h = h(x) 

Fh = J g ■ hdx. (46) 

This gives us access to the representation (1251) of its Frechet derivative 
considered as a nonlinear functional FVcrijk = F(a ) of a. Using this repre¬ 
sentation, a partition into dyadic annuli, and Holder’s estimate (recall (1HT]) ) 
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we obtain 


\ dF \\ 
da 9 

< ||(|Vu| + |Vtij-feDIV0j + e*| || g ,| x |<2r 

oo 

+ £ll(|V«| + |V^|)|V</, i + 

e i\^q,2 n r<\x\<2 n + 1 r 

n —1 

(|| ||p i | iC |<2r “1“ || ||p,|cc|<2r) || ^ || 2,\x\<2r 

oo 

+ £(liv»||* n r<|rr|<2 n + 1 r II ||p,2 n r<|a;|<2 n + 1 r )l|v& + ^i\\2,2 n r<\x\<2 n + 1 r- 

n=l 

(47) 


In view of (j32j) applied to the a-harmonic function u{x) — x t + (pi{x), cf. (|3|), 
we obtain for all radii p > r using Caccioppoli’s inequality and 

/ | V0j + ei\ 2 dx < p d lim inf 4 - | V0, + ej| 2 da; 

J |a;|<p ^1°° i|x|<_R 


< lim inf 


R t°° R 2 J\x\<2R 


k + Xi 


(0i + Xi) 


p|<2i? 


dx < p d . 


Hence (pT71) turns into 


9F \\ 
da 9 


< rH\\Vv\\ p +\\Vv jk \\ p ) 


OO 

“t - ^ 2 (|| ||p,|:c|>2 n r “t - || ||p,|a:|>2 n r‘) • 

n=l 


(48) 

(49) 


It thus remains to estimate the auxiliary functions v and Vjk■ The estimate 
of the terms in line (HSjl is easy: By (H51) and Calderon-Zygmund for (124)) 

we obtain ||Vu|| p < ||g|| p < r p . By (141J) we have Calderon-Zygmund 

with exponent p for the equation (125)1 . so that || Vu^Hp < ||Vu|| p < r ~ . 
In order to control the terms in line (1491) . we shall establish the following 
estimates for n E N 

II V^||p,|x|> 2 "r < (2 n )~ d+ ir- E ^ d , (50) 

W^VjkWp,\x\> 2 ^r < n{2 n )~ d+ Pr~ R ^ d . (51) 

We note that since p > 2, these estimates imply that the sum over n in f!49p 
converges and gives (IMP . 
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The estimate (15U]I for the solution v of the constant coefficient equation (PH) 
is classical: We already argued that || Vu|| p < r ~ , by the estimate on the 
support of g in (l45li we have that v is harmonic in {|x| > r} and that it has 
vanishing flux f,,x ■ S7v = 0. It thus decays as |Vv(x)| < \x\~ d r d ~p || Vu|| p 
for |x| > 2r, which in particular yields fl50p . We now turn to (15ip and to this 
purpose rewrite the equation (1251) for Vjk as 

— V • a*X7vjk — V ■ g 

with the r. h. s. g := —a*(djVek — dkvej). We already argued that ||g|| p < 
r p and ( I50|) translates into 


p,\x\>2 n r 


< 


(2 r 


v-d+- - 
) Pf 


-2=1 C 

P 


(52) 


In order to proceed, we split g into {g m }m =where g 0 is supported in 
{|x| < 2 r} and for m > 1 g m is supported in { 2 m r < |x| < 2 m+1 r}, so that 
([52]) translates into 

\\9m\\ p < (2 ™y d+ ir-^ d . (53) 

This entails a splitting of Vjk into {u m } m =o,i ,- 5 where v rn is the Lax-Milgram 
solution of 

— V • a*Vh m = V • g m . (54) 

We will now argue that 


||Vfi m || p , w > 2 n r <min{(2")- d+ p,(2 m ) 


' p 


}r p , 


(55) 


which implies the estimate (T5Tj) by the triangle inequality || yvjk\\ p ,\ x \> 2 n r < 
J^m=0 II We note that ([55]) together with our Calderon-Zygmund 

estimate (jUD applied to (1541) yields ||Vu m || p < (2 m ) d+ pr P p d . In order to 
establish (1551) . it thus remains to show 

l|V6m||p,[x[>2«r ^ (2 n )~ d+ pr~ E ^~ d for m < n. (56) 

p 

We argue in favor of (1561) by duality and thus consider an arbitrary h e 
supported in {|x| > 2 n r} and denote by w the corresponding Lax-Milgram 
solution of (FIUp . By integration by parts, we deduce from (FlUj) and (1541) that 
f h ■ 'Vv m dx — f g m ■ Vwdx. By the support condition on g m this yields 


Jh.Vr wlx 


— Ilfl'mllpH Vic|| —JL_ ; |ir|<2 m + 1 r • 
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By the support assumption on h we have that w is a-harmonic in (|x| < 2 n r}. 
Since m < n, we may use (02]) applied to w in form of 


l|V“ll*,W< 2 ”+ir 5 (2”- m )-' + 5||Vtt>|| AiM£2 


We combine this with (I4TT) in form of ||Vrc||j_ < ||/z||_e_, and with (]53|) . to 
obtain 


h ■ Vv m dx 


< 


(2" 


d+i -E=ld 

) Pf V 


P-1 


which gives (156]) . 


In the case of a functional of the scalar potential of the form F(a ) = 
we claim that the Frechet derivative of F is again controlled in the sense of 
(|34p . The proof is mostly analogous to the previous one; we again rewrite F 
as in (1461) with some g satisfying (1451) . Starting from the representation (1281) . 
one derives an analogue of estimate (1471) reading 


dF 
da q 

~ || Vu||p,|a;|<2r ||T Cj 112,|a;|<2r- 

oo 

+ II^llp,2"r<|a;|<2"+ 1 r|| V0i + 61| 2,2 n r<|a;|<2 rl + 1 r • 
n =1 


The second factors on the right in this estimate coincide with the ones in the 
case F(a) = FVa^k] therefore, we get the following analogue to estimate 

m- 


dF 
da q 


< r*||VU|| p + £(2V)*||W|U|>^. 

n= 1 


The equation (1271) for v has the structure of the equation (1541) with m — 0 
(including the estimate ||c/||j, < r p~ and the inclusion supp g C {|x| < 2r}, 
cf. (051) ); therefore, the decay property (156|) carries over to our v. This 
establishes the estimate 


dF 


<r -*^<1 

\ nr 2 p 


also in the case F(a) = F’Vfa. 


□ 
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3.4 Sufficient conditions for the mean value property 
in terms of linear functionals of the corrector 


Proof of Lemma\ff In order to show that (j36j) and (T33l) imply (l32]h we only 
need to show the existence of functionals F\ , • •' ,Fn such that (T36|) and (133|) 
imply 


1 

R 



K0.O-) 



•C 1 for all R> r. 


(57) 


By Proposition [21 the estimate (122]) follows from (15711 . 

Let us now give the argument for fl57jh First, it is clearly enough to show 
that for any 0 < S 1, there exists N < 5~ d functionals iy, ■ ■ • , Fv on 
vector fields which are bounded in the sense of 


P-1 


\F n h\ <r d ( / \h\-^dx) ? 

J |a;|<l 


(58) 


and such that for any dyadic p > 1 


P ' J \x\<p 




hi <p 


\ Z \ 2 

(0,<r) efo) <5+ sup max \F ntR V{<j>, a 


R>p dyadic n ~ L”' 


By dyadic iteration, it is enough to show for any dyadic p > 1 


.. K^ 0- ) ■ 

P'd\ x \<p J\x\< P 

< max \F np V (0, cr)| 

n= l,--- ,N 


(0, a)\ 2 dx 


Kb 


\x\<2 P 


lOM-f ( 0 , 0 -)| dx) +1 


\x\<2p 


(59) 


(60) 


Indeed, abbreviating D m := -, <2m |(0, or) - f 2m (0,a)| 2 dx) 2 , the es¬ 


timate PH) may be rewritten as (using a slight readjustment of 5) 


D m < C 0 max |F n 2 mV(0, a)| + S[ D m+1 + 1), 

71=1,— ,JV V / 

which may be iterated to 

Dm < 1 Co max max |F n 2 MV(0, cr)| 

1 — 0 ,m+mo ,N 

+ + 5 mo+1 D m+mo+1 . 

l — o 
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By our sublinearity assumption on the corrector (fool) (which may be rewritten 
as lirn motoo D mo = 0), this yields (E2D- 

We now turn to the argument for (l60j) . By Caccioppoli’s estimate on (J4]) we 
have 


M <|p 


\\7(f>i\ 2 dx 


1 < I 

rsj 


P ' J\x\<2p 


(<j>i 


\x\<2p 


< t>i) Z dx) 2 + p d/2 , 


and thus in particular for the flux g* = a(V</g + e*) 

1 


' W <fp 


|g;| dx 


< 


P ' J\x\<2p 


\x\<2p 


•i ) dx) " + p' 


d/2 


Caccioppoli’s estimate on ((9]) gives 




1 

2 


< 

r 


< 


( / ( a ijk 

V d|x|<|p 

/ a ijk ) 2 dx) 2 + ( 

7 M<§ p ' V 

[ \qi\ 2 dx)' 

■J\x\<Ip 7 

( / { a ijk ~ ' 

'J\x\<2p 

/ a ijk ) 2 dx) 2 + ( 

l\x\<2 p 7 v 

[ |g*| 2 dx) 2 

J\x\<lp 7 


The last three estimates combine to 



Hence for (1601) is enough to show 


\x\<2p 


(0, a) dx) + p' 


d/2 


1 

P 


(/ | (</>,u)| 2 dx ) 2 

V|a;|<p d|x|<p 7 

< max |-F np V((/>, cr)| + sfi- |V(0,cr)| 2 dx 
n=1 ’~’ N ’ V|x|<p 


(61) 


This statement is not just true for (</>, a) — jj x ^ <p ((p, a), but for any function 
( of vanishing spatial average on {|x| < p}: By rescaling, it is sufficient to 
show the estimate on the unit ball {|x| < 1}. It is more convenient to see it 
when the unit ball {|x| < 1} is replaced by the unit square (0, l) d : 



< 


max \F n X7(\ 2 + S 2 [ |VC| 2 dx. 

n=l,-, N J (01)d 


(62) 
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Indeed, dividing (0, l) d into N = 6 d (suppose that 5 1 is an integer) sub¬ 
cubes {Q n }n=i, -,N °f side length 5 and setting F n V( := J ^ (dx (recall that 
/( oi) d ^d x = 0 so that F n is indeed a function of V((), (1D5D follows from 
using Poincare’s estimate on each Q n in form of L ( 2 dx — | QuH^q (dx) 2 < 
5 2 | VC | 2 dx and then summing up. We note that by Poincare’s estimate 

on (0, l) rf , the F n have the desired boundedness property (1551) . at first on 
gradient fields V( 

lF n V(j<S~ d ( [ |C| ^dx)^ <5~ d ( [ |VC| ^dx)^, 

J( 0,l) d 4(0,l) d 

and then on any vector held h by extension a la Hahn-Banach. □ 


4 Proof of Main Result 


Proof of Theorem\l\ 

Proof of Assertion i). 

Consider the functionals {-F n }n=i,- ,n and their rescalings {F n , R } n= i,- ,n-,r dyadic 
constructed in Lemma 0J Let us abbreviate F ni ^V(0, cr) as F n R . We would 
like to apply concentration of measure to these functionals. 

The main difficulty that we need to overcome is that our sensitivity estimate 
(134)) in Lemma [3] for the quantity F Utr is based on the assumption that the 
mean-value property (1321) holds for a-harmonic functions down to scale r. 
By Lemma [4] this assumption may be reduced to the smallness assumption 
(1551) for our functionals F m)R on scales R > r, so that Lemma [5] becomes 
applicable under the assumption (1551) : Let q be related to fd through (1511) 
and let p be related to q through (151]) . By the smallness assumption on fd in 
our theorem (cf. (1151) ). we deduce that (1501) holds. By scaling, our functionals 
F n>r satisfy the estimate (1551) up to a universal constant factor. Furthermore, 
by ergodicity the property (1551) holds for (-)-almost every coefficient held a 
(regarding a, this result has been shown in [ID] Lemma 1 ]; for </>, it is classical 
but may also be found in [ID]). Thus, the estimate (1541) holds for F nr under 
the assumption (156]) . i.e. there exists a constant Co only depending on d, A, 
and /3, such that for any n — 1, ■ ■ • , N and any radius r the implication 


sup 

m,R>r dyadic 


Fm,R < 


1 

Co 


,dF, 


da 


n ' r m 2 < 

q 


1 

yd 


(63) 


holds for (-)-a.e. coefficient held a. 

To apply concentration of measure in the form of Proposition [Tj to some 
functional F, we however need an unconditional bound on the Malliavin 
derivative (cf. (ITS]) respectively (155]) ). 
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Therefore we first introduce a new random variable whose derivative vanishes 
whenever the smallness condition in (Rn?j) is violated: Consider the auxiliary 
random variable 

F r '■= min{ sup \F m<R \, —}, (64) 

m,R>r Lq 


where the sup runs over all dyadic radii R = 2 k r, k G No- By the usual 
differentiation rules applied to the Frechet derivative ^ in the norm || • || q , 
we obtain 


dF r 

da 


< /( sup F, 

m,R>r 


m.R 


1 . 

< 7 -) sup 

C'O m,R>r 


dF m , R 


and thus by (163|) 

,,dF r 2 < 1 

By Lemma [H we may apply concentration of measure in form of (1201) to the 
random variable crF 2 F r (where c is some small universal constant). This 
yields 

(I(\F r - (F r )| > M)) < exp(--r /3 M 2 ) for all M > 0, (65) 

o 

so that it remains to control the expectation (F r ). 


Because of (jHJ) and the definition of F m>R , it follows from qualitative ergodicity 
of (•) and Birkhoff’s ergodic theorem that lim^oo F m R = 0 almost surely, so 
that by dominated convergence \im r ^ 00 (F r ) = 0. ffence there exists a finite 
radius ro which is minimal with the property 


(Fr) 



for all r > tq. 


( 66 ) 


ffence using (jf?5|) for M = jf- we get in view of the definition f)64p of F r 


(/( sup \F m<R 

m, R>r 


> 


l 

2 (To 


)) ^ exp( 



for all r > tq. 


(67) 


On the basis of (l67lh we now get a quantitative estimate on r 0 . To this 
purpose we now consider the auxiliary variable 

Fn,r ■ f]( sup F m R )F nr j (68) 

m,R>r 

where again the sup runs over all dyadic radii R = 2 k r, k G No, and where 
the cut-off function r] = 7](F) is given by 

rj(F) = max{min{2C' 0 (—- F), 1}, 0}. (69) 

Co 
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The advantage of the auxiliary variable (1551) over (15-ID is that we control 
its expectation: Since the stationary V(</>, a) has vanishing expectation, cf. 
(0), and by the linearity of F n>r in V(0, a) we have (F Htr ) = 0 and thus 
{F n , r ) = {(v — 1 )F ntr ) so that by construction of r] 

K-Fn,r>| < (/( sup \F m , R \ > -^)\F n , r \). 

m,R>r ^^0 


Since the stationary V (0, a) has bounded second moments, cf. flSJ), and by the 
boundedness property of F n r in V(</>, a) we obtain from the Cauchy-Schwarz 
inequality 

(F n ,r ) 2 ^ U( sup | F mtR \ > ^r)), 

m,R>r ^^0 

which in view of (157D improves to 

I (K,r) | < exp(-^r^) for any r > r 0 . (70) 

By differentiation rules for the Frechet derivative F in the norm || • || g we 
obtain for the auxiliary random variable F n:1 . 


,dF n 


da 


sup j 

m,R>r 


< /(sup \F m>R \ < —)(2C 0 \F n 

m,R>r U-O 

^ nr/ I Tp I ^ I \ II dF m R 

< 3/( sup \F mtR \<—) sup ||——||„ 

U (} m,R>r Cla 


, dF m ,R .. ,, dF n , r n \ 

da q da q ' 


m.R>r 


and thus by (153]) 


dF n r 2 < J_ 
da 1,9 ~ rF 


and hence by concentration of measure in form of (12U|) (applied to cr^Fp n r 
by means of Lemma HI c being a small universal constant) 


(I(\F n ,r - (F n , r )I > M)) < exp(-Ir^M 2 ) for all M > 0. 
Together with (170D this yields 


(J(|F n>r | > M)) < exp(—-^r^M 2 ) for M exp(—-^r^) and r > r 0 . 

Cs G 

By definition (J55|) we have {I{\F n)T .\ > M)) < (/(sup mR > r \F m)R \ > FF)) + 
(I(\F ntr \ > M)) so that by (157)) the above upgrades to 


(/(|F n , r |>M))<exp(-^M 2 ) 


for 1 > M » exp( 



and r > tq . 
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Since exp(— ^r 13 ) < 1 for all r, the above holds without the lower restric¬ 

tion on M: 

(I(\F n ^ r \ > M)) < exp(—for M < 1 and r > r 0 . (71) 

o 

Using this estimate with r replaced by R and summing over the finite index 
set n — 1, • • ■ , N and all dyadic R > r we obtain 

(/( sup |F m ,*| > M)> < exp(-ir^M 2 ) 

m,R>r ^ 

for M < 1 and r > r 0 (72) 

and thus in particular for the auxiliary random variable (164)) 

(I(F r > M )) < exp(——r /3 M 2 ) for all M > 0 and r > ro, 

o 

where the upper bound on M is immaterial since F r < ^ < 1. Using (F r ) = 
J 0 °°(/(-F r > M))dM, this yields the following quantification of lim r ^ 00 (F r ) = 
0: 

/»oo 

(F r ) < / exp (—r 13 M 2 )dM < r _ 2 for all r > r 0 . 

Jo 

Since r 0 was minimal in (15^1) and since (-F r ) depends continuously on r, this 
yields the desired 

ro < I- (73) 

It remains to argue why CD, which together with fl73|) may be rephrased as 

(J(|F nir | > M )) < exp(—-^r^M 2 ) for M < 1 and r 1, 

yields ([13]) . It just suffices to include the given functional F from (1X3]) into 
the list of finitely many functionals Fi, ■ ■ ■ , F,v, say, as the last functional 
Fn = F, and then to specify the above to n = N. We note that for q related 
to /3 through (l2Tf and p related to q through (l3Tj) one has i.e. 

m entails (j29j) . Note that by adjusting the constants, ( 1151) is trivial for 
r < 1, so that we obtain (fT7>l) over the whole range r > 0. 

Proof of Assertion ii). 

The arguments in this section require /3 < 2, which in view of our assumption 
[3 -C 1 is no restriction. Let r* denote the minimal dyadic radius with the 
property (TT6l) - we know but don’t have to use that it is finite by quantitative 
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ergoclicity. In order to establish (TT7T) . it is enough to show for a given dyadic 
r 0 > 1 that 

(/(r* > r 0 )) < exp(-^rj). (74) 

It will be convenient to replace balls by cubes. Moreover, all radii or rather 
side length are dyadic. By definition of r* as the smallest radius with fflfiH . 
the event r* > ro means that there exists a radius R > r o with 


1 

Jr 


(~R,R) d 





(75) 


where 


f(z ) := log(e + log z). 


In the sequel, the intermediate (dyadic) radius r 1 e [r 0 , R] with 


T\ 




< r A\P 

K R 


/( 


R, 

'r 0 ' 


(76) 


will play a role. Note that we use here f3 > 0 and that f(z) grows sub- 
algebraically. For the 1. h. s. of (T75|) we note 


(~R,R) d 


K&o-) - 


(~R,R) d 


(0, a)\ 2 dx 


(~R,R) d 


|(0,cr) - (0,o-) ri | 2 da:+ 


re[2ri ,R] J (~^R) d 
r dyadic 


I 2 d.X, 


(77) 


where (0, cr) r denotes the L 2 ((—i?, i?) d )-orthogonal projection of (0, ex) onto 
the space of functions that are piecewise constant on the (—) d dyadic sub¬ 
cubes Q of “level r” (that is, of side length 2r) of the cube (-R, R) d ■ In other 
words, on such a sub-cube Q, (0, a) r = $q(4>i cr)dx. With this language, we 
may rewrite the first r. h. s. term of (177)1 as 

|(0W) - (0,^)rj 2 dx = (J±) d ^ j- l(^> a ) (<^)| 2 dz, 


(~R,R) d - q level ri • 

so that by Poincare’s estimate on each of the cubes Q we obtain 


{-R,R) d 


1(0, a) - (0, cr) ri | 2 c7r < rl -f | V(0, <r)\ 2 dx, 

J{-R,R) d 
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and then by Caccioppoli’s estimate based on (3} & (EJ), cf. (RTTll . 


(~R,R) d 

<(-) 2 


|(0,a) - (0, <y) ri \ 2 dx 

l(0w) - / 

(—2R,2R) d J(-2R;2R) d 


(0, a)\ 2 dx + r 2 . 


As a consequence of LemmaHl there exist N ~ 1 linear functionals {i^n}n=i,- ,n 
whose rescaled versions F n/r satisfy the boundedness property (TT4|) such that 
for any r > 2 R we have the implication 


sup max (i 7 l„ )? .V(0, cr)) 2 <C 1 

r>2R dyadic n= b"'A 


1 

1$ 


1(0; a ) ~ 

'(-2R,2R) d 

From the two last statements we gather 


(~2R,2R) d 


((p,a)\ 2 dx < 1. 


sup max (F n>J .V(0, a)) 2 -C 1 

r>2R dyadic n=1 -'" A 

|(0,a) - (0,a) n | 2 cte < (^) 2 . 
O J(-R,R)d R 

In view of (17(1 this can be rewritten as 


V r > 2i? dyadic, n = 1, 
1 


,iV (F n , r V(0,a)) 2 <l 


A 2 


(~R,R) d 


1(0, a) - (0, cr) ri | 2 dx < ^(^) /3 /(^); ( 78 ) 


provided that we adjust the definition (1751) of rq appropriately (to obtain the 
estimate < in (1781) in place of just <). We now turn to the second r. h. 
s. term in dZZD, which in view of the definition of (0, cr) r we may estimate as 
follows 


Y f I(0W)§ - {<f>,cr) r \ 2 dx 

re[2ri,-R] J(~ R ’R) d 
r dyadic 


< 


£ 

re[2ri,_R] 
r dyadic 


max max 

Q level r Q'cQ level | 




Hence if for any of the (^) d dyadic sub-cubes Q of (—R,R) d of level r we 
introduce the N = 2 d linear functionals Fq n as an extension of 

Fq.t^C ■= -(/ C dx - / C dx), 
r JQ'n JQ 
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where {Q' n } n=1 ... 2 d is an enumeration of the sub-cubes of level | of Q, and 
which satisfy the desired boundedness property (ITT]) restricted to gradient 
fields (which is no issue because of Hahn-Banach extension) and translated 
(which will be no issue because of stationarity), that is, 


l*Q,nVC|<(/ |VC| ^dx) d ™ : 

JQ 


we have 


Xi T (<#•. <d§ - (<t>,v)r\ 2 dx 

r dyadic 

< (4) 2 max (F Qin V(0,or)) 5 

z ' R Q level r n =l,— ,2 d 
rG[2r u R] 
r dyadic 


(79) 


From this we learn, since for the auxiliary function g(z) := log 2 {z + e), the 
dyadic sum X]re[ 2 n r] d(^) i s universally bounded, 


V r G [2r 1 , R] dyadic, Q level r, n = 1, • • • , 2 d 

(F«,„vW^)) 2 «(-) 2 s(-)(^)' , /(-) 

r r R r 0 


1 


V 

re[2ri,R] R~ R ^ R ) d 


1(0,a)r - (0,cr) r | 2 da; < 


r dyadic 


In view of (1771) . the combination of this with (178]) yields 
V r e [2r 1 , R] dyadic, Q level r, n — 1, ■ ■ ■ , 2 d 
(F Q ,nV(M? 

r r R r 0 

and V r > 2R dyadic, n — 1, • • • , iV (F„ r V(0, a)) 2 1 


1 


{-R,R) d 


1(0, Or) - 


(0,a)| 2 oh < (^) /3 /(—)• 


/?. 


(-R,R) d 


R' 


R 0 


(80) 


Equipped with this deterministic argument, we now may proceed to the 
stochastic part: In the event of r* > r 0 , there exists a dyadic R > r 0 such 
that (1751) holds, so that we learn from (180|) that there exists 

• a (dyadic) r e [2ri, R], a sub-cube Q of (—R,R) d of level r, and an 
index n= 1, ■ ■ ■ , 2 d such that (F Q>n V(0, a)) 2 > (f) 2 ^(f)(f ) /3 /(^)- In 
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view of the boundedness condition (1791) and stationarity, we may apply 
(113]) with F replaced by Fq jH and M 2 replaced by (y) 2 #^)^)^/^)- 
This M is admissible in the sense of M < 1 because by (176]) we have 
(*) 2 9(7)(i£ (^) 2 (3)' 5 /(5) ~ 1- Hence the probability of 
each single of this events is estimated as follows 

</((f q ,„v«.,o )) 2 > i<f 

< exp (-^(-) 2 -V-)/(-)>f)- 

C r rro 

Since g(z) decays sub-algebraically in z and since j3 < 2, this yields the 
simpler form 

{l((F Q , n V(M? > 5(f) 2 9(f)(|)"/(^))> 
iS exp(-i(^)‘-?/(^) ro s ). 

• or a (dyadic) r > 2 R and an index n — 1, • • • , N for which the estimate 
(-Pn,rV(0, a)) 2 > 1 holds. By the boundedness property of F n . r , each 
single of these events is estimated as 

<f((f«,rV (0, cr)) 2 > ^)> < exp(-^r /3 ). 

Taking the number (y) d of sub-cubes Q into account and recalling N < 1, 
this implies 

(J(r„ > r 0 )) 

s E ( E (7)' , “p(-^7) 1 - s /(^)^)+ £ exp(-i r 9 )). 

R>ro rG[2ro,i?] ^ r>2R 

R dyadic r dyadic r dyadic 

(81) 

Again, since 1 — ^ > 0, we have the calculus estimate 

rE[2ro,.R] 
r dyadic 

<exp(-yl) y (^exp(-IrUo g (^)) 

re[2r 0 ,fi] 
r dyadic 

< exp(—A) for A>1. 
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Applying this to the first sum over r in (I8ip and A = which sat¬ 

isfies A > 1 for r 0 > 1, and using the estimate Yl r >2R- r dyadic exp(-^r /3 ) < 
exp(— -^R 13 ) (which holds provided that R > tq > 1 ) for the second sum over 
r, we obtain 


(J(r*>r 0 )) < 



R dyadic 


for r 0 ^ 1. 


Thanks to /3 > 0, we have exp(—< exp(—^/(^)^o)> so that the second 
summand is dominated by the first one: 


(J(r*>r 0 )) < 



R dyadic 


Now we see the reason for the choice of f(z) = log(e + logz) for which 
/(2 m ) > -|- log(m + 1)) and thus 

oo 00 A A 

^exp(-A/(2 m )) < ^exp(--(l + log(m + l))) < exp(--) for A » 1. 

m =0 m =0 

With ^rg playing the role of A this yields ( 174 ]) . Note that the condition 
r 0 1 is immaterial after adjusting the constants, as the 1. h. s. of ( 1741 ) is 
bounded by 1. □ 
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